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Wilson's area law in QCD is critically examined. It is shown that the expectation value of 
the Wilson loop integral exp(J iA^dx 11 ) in the strong coupling limit vanishes when we employ the 
conjugate Wilson action which has a proper QED action in the continuum limit. The finite value of 
Wilson loop with the Wilson action is due to the result of the artifact. The fact that his area law 
is obtained even for QED simply indicates that the area law is unphysical. 
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I. INTRODUCTION 

Since Wilson pj proposed Wilson loop in evaluating 
the gauge field theories on the lattice, it has been com- 
monly accepted that Wilson's area law is the basic con- 
ceptual ingredients to understand a confining mechanism 
in QCD. Wilson's formulation of the lattice version of 
QCD is interesting, and in addition it can be used for 
evaluating the lattice QCD in terms of numerical calcu- 
lations. Indeed, there are many calculations of the lattice 
QCD simulations which confirm the confining mechanism 
of Wilson's criteria 

In his paper, Wilson presented an area law which can 
suggest a confining potential of the following type Q] 



V(R) 



ln.g 2 



R 



(1) 



where a denotes a lattice constant. Therefore, it means 
that the QCD confining potential is a linear rising one. 
Further, this shape of the potential is well reproduced by 
the lattice QCD simulations |3|. 

However, one may have somewhat an uneasy feeling on 
eq.(l). That is, eq.(l) does not have a proper dimension 
since the lattice constant (units) cannot play a role of a 
physical quantity. In fact, if we take the continuum limit 
of eq.(l), then the right hand side goes to infinity, and 
therefore it becomes unphysical. Here, we should explain 
the continuum limit itself in order to avoid any confu- 
sions. This is closely related to the way one solves field 
theory. When we wish to solve field theory models, the 
space and time are of course from — oo to oo. However, 
we cannot solve the field theory models in this space and 
time, and therefore normally we put the theory into a box 
with its length L. In this case, we can solve field theory 
models within the box and after calculations we should 
make the length L much larger than any other scales of 
the models, which is called thermodynamic limit. If this 
continuum field theory model can be solved in some way 
or the other, then we can obtain physical observables af- 
ter the thermodynamic limit is taken. 



Sometimes people want to solve the field theory by 
cutting the space and time into a lattice. Wilson's way of 
solving QCD is just this lattice field theory. In this case, 
we should cut the space and time by the site number N, 
and therefore the lattice constant a becomes 



L 

N 



(2) 



where it is to be noted that the site number N is not a 
particle number of field theory model. In this case, it is 
clear that the lattice constant is not a physical quantity 
since it depends as to how one wishes to cut the space 
and time. In fact, one sees that the potential [eq.(l)] 
indeed diverges when one wants to cut the space as small 
as possible to simulate the continuum space, keeping the 
box length L finite. Therefore, it is clear that, physically, 
eq.(l) does not make sense. 

Here, it should be important to note that the lattice 
constant a in solid state physics has a completely differ- 
ent meaning from the continuum field theory since the 
constant a of the lattice is fixed and finite. In this case, 
therefore, there is no ambiguity of cutting the space, and 
thus the thermodynamic limit is the only concern in the 
calculation of physical observables. 

In this paper, we show that the finite value of the Wil- 
son loop is an artifact due to the Wilson action since it 
should inevitably pick up unphysical contributions to the 
loop integral. A similar work has been done by Grady 
who shows that the action different from the Wilson 
action gives no area law in SU(2) lattice gauge theory. 
There are some works which critically examine lattice 
calculations 0,0; but we do not go into details here. 



II. AMBIGUITIES IN AVERAGING 
PROCEDURE 

Now, we should discuss the ambiguities in evaluation of 
the expectation values. The discussion here will be used 
in the later sections to understand the way of evaluating 
the average of a Wilson loop. As examples, we employ 
the distribution function p± (A) as 
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p±(A)=exp (he ±lA )=J2 



+ikA 



k=Q 
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(3) 
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where h is a constant. This type of the distribution func- 
tion is often used in mathematics in order to pick up 
some integer from certain complicated functions. 

Now, we make an expectation value of A and obtain 



oo , k f x 



where X should be eventually set to infinity. Eq.(4) 
shows that there is no finite expectation value of A with 
the distribution function of p±{A) in eq.(3). 

However, if we make expectation values of e mA where 
n is a positive integer, we obtain 



(e mA ) + 



(e mA ). 
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P-(A)dA = — 
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(6) 

The expectation value of (e mj4 }_ survives even though 
the expectation value of (e mA )+ vanishes. The last re- 
sult [eq.(6)] shows just the same trick as Wilson used in 
obtaining a finite expectation value of the loop integral, 
and one notices that the procedure is just to pick up the 
integer value n in the function of e mA . 



ILL 



WILSON'S AREA LAW 



A dimensionless form of the field strength F njiU is defined 



as 



{A)±=J2ji Ae ±lkA dA/ P± (A)dA = (4) vector field A, 



fn[iv — ^ Q^n^iv — B n ^ -\- Bn+pi,u B n -^{j^ B nu . (10) 

Wilson imposed the periodic boundary conditions on the 



where N is related to the box length L by L = Na. 
In the strong coupling limit, eq.(7) can be expressed 

as 



nn/ dB ^ e e 
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In the same way as eq.(6), eq.(ll) has a finite contri- 
bution only when the exponent in the integrand is zero. 
Therefore, the nonzero terms in the sum are those for 
which satisfy 
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(12) 



Now, we come to the area law in Wilson's paper. Most 
of the lattice formulations of the field theory are well 
written in his paper, and QCD on a lattice itself is quite 
interesting. In particular, he introduced Wilson's loop 
and evaluated it in terms of Euclidean path integral for- 
mulation. Here, we repeat the discussion in his paper 
using the same notations. The expectation value of the 
Wilson loop is written as 



(exp^O )+=^ i nn /" 



dB 



nifi 
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xexp ^B lltl + l^ e ^ J 
where the partition function Z + is defined as 

z + = nn f dB ^p f^E^) 



(7) 



(8) 



Now, if we specify the path P which contains K pla- 
quette, then we should find the number of fi k . Kk a k with 
k = K. In this case, we see that eq.(12) is indeed satis- 
fied. Therefore, the k is the number of the plaquette that 
are surrounded by the path P. If this area is denoted by 
A, then wc find that k = K = Ar. 

Therefore, it is easy to find 



exp i^2,B ntJl )_ 



(5 2 ) 



-A/a 2 



(13) 



where A should be described by some physical quantity 
like A = RT with T time distance in Euclidean space. 

From eq.(13), one obtains eq.(l) using the method of 
transfer matrix which relates the Wilson loop to the po- 
tential 0,0. However, eq.(13) has a physically improper 
expression. That is, the area A is described in terms of 
the physical quantity while a 2 is not a physical quantity 
since it should be eventually put to zero. Therefore, it 
is by now clear that the right hand side cannot survive 
when one takes the continuum limit. 



P denotes a path on the lattice. Here, the dimensionless 
quantity B nfl is related to the discretized vector potential 



A n ^i as 



B n ^ = agA n 



(9) 



IV. 



A TRICK IN WILSON'S CALCULATION 



Why did he obtain such an unphysical result ? The an- 
swer is simple. The problem is related to the ambiguity 
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in the definition of the partition function (8). Wilson's 
area law is obtained only when the Wilson action is em- 
ployed. All other actions including Grady's action Q 
cannot reproduce the area law. We will show it below by 
a simple action with a negative sign in the exponential. 

In eq.(8) the gauge field action is expressed as 
3+ = "Tg 1 Eniu- e^ n "" instead of the conventional form 
— j a 4 X^rmi/ ^np,v For small a, this term can be expanded 
as 

s + = h £ ei/ " M " « h D 1 + ^ \&» ■ ■ ■)■ 

The first term is irrelevant and the linear term in f niiv 
gives no contribution in the calculation. The quadratic 
term corresponds to — j^^n/iu^n^ which is the con- 
ventional expression of the gauge field action. Here, one 
can easily notice that it is also possible to write the action 
for the gauge field with a minus sign in the exponential 
as S- = J2 nf iv e -4 ^""" . Because S- also reduces to 

5_ ~ —\a 4 Ylnuu Fn)iu f° r small a. Therefore there is an 
ambiguity in the definition of Z . 

The partition function defined with S- is written as 

In the strong coupling limit, the expectation value of the 
Wilson loop with Z- becomes 

*im x; - e 

m M " £i7Ticri Ik^k^k 




As in eq.(ll), the non- vanishing terms in eq.(16) are 
those for which satisfy 

^ Brit* ~ //uriffi fe k K k o-k = 0- (17) 

p 

However, for the same path P as in eq.(ll), there is no 
term in eq.(16) that satisfies the above condition. Any 
set of — fg n <j'm eq.(16) cannot produce a loop in the oppo- 
site direction to the path P, and therefore cannot cancel 
J2pBntt- The situation is the same as eq.(5). Thus the 
average of the Wilson loop vanishes to zero 

<expfi^B n ^)_=0. (18) 



Eqs.(13) and (18) show that for a given path P 
{exp(iJ2 P B nf _ l ))+ and (exp(i J2p B n»))- behaves differ- 
ently. This is because it is impossible to satisfy both 
eqs.(12) and (17) simultaneously for a given path P. 

Here, we should make a comment on the choice of the 
action Wilson employed. There is some argument that 
the action should be Hermitian, and therefore one should 
take the action of (5+ + 5_)/2 which is indeed Hermitian. 
However, this modification of the action is a matter of no 
significance since the important part of the action is of 
course —j^J^n^-^n^ f° r small a. Therefore, if the 
Hermiticity of the action is an important factor in the 
evaluation of physical observables, then it indicates that 
the action of eq.(14) as well as the Hermitian action is 
physically not a proper one. 

Further, if one calculates the average using the con- 
ventional action of the gauge field, —ja 4 J2 n ^ v F^ v , 
then one obtains a different result which has no area 
law. In addition, if one employs an action, for example, 
S' = g|s Yln/xv e l2 ^ n ^ , which has a proper continuum 
limit, then one obtains no contribution to the Wilson 
loop integral. 

From the above discussion, it is concluded that the 
Wilson loop integral [eq.(13)] is accidentally finite as an 
artifact of the averaging procedure. This also indicates 
that the potential V[R) ~ obtained from eq.(13) 

has nothing to do with physics, which confirms our point 
in the first section. 



V. WHAT IS THEN CONFINEMENT 
MECHANISM 

We have proved that Wilson's area law has nothing to 
do with the confinement in QCD. The basic mistake of 
Wilson's calculation is due to his action which inevitably 
picks up unphysical contributions even though Wilson's 
action can be reduced to a proper gauge field action in 
the continuum limit. However, higher order terms of his 
action contribute nonperturbatively to the Wilson loop 
integral and generate a fictitious area law as a result. 
Therefore, the "confining potential" which has been be- 
lieved to be confirmed by the numerical calculations of 
the lattice QCD should be reexamined. 

From the dynamical point of view, it is obvious that 
gluon dynamics without fermions cannot give the con- 
finement mechanism since gluons produce interactions 
between quarks which should be always Coulomb like po- 
tential. Thus, one has to consider quark and anti-quark 
pairs in the intermediate states since they have masses 
and thus can induce a new scale. Therefore, physics of 
the confinement should be obtained only after one solved 
the QCD dynamics with quarks and gluons together. In 
this respect, when one carries out the lattice simulations 
of QCD, the quenched approximation is clearly meaning- 
less. One should take into account the fermion degrees 
of freedom properly. 

Now, what should be a new mechanism of confining 
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quarks ? The confinement itself must be closely related 
to the gauge invariance of physical observables. The 
color non-singlet objects should be gauge dependent, and 
therefore they cannot be physically observed. The quark 
confinement should be due to the gauge non-invariancc 
of the color charge and, because of it, the dynamical con- 
finement should occur accordingly. In this respect, the 
fact that Wilson's area law is proved for the U(l) gauge 
theory as well clearly shows that his calculated result 
should not have any relevance to the confinement mech- 
anism in QCD @. 



In this sense, it is most probable that the essence of 
the confinement may well be realized by the MIT bag 
model H3. 
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